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Abstract 

We present a phenomenological application of the joint resummation formalism to elec- 
troweak annihilation processes at measured boson momentum Qt- This formalism simul- 
taneously resums at next-to-leading logarithmic accuracy large threshold and recoil correc- 
tions to partonic scattering. We invert the impact parameter transform using a previously 
described analytic continuation procedure. This leads to a well-defined, resummed pertur- 
bative cross section for all nonzero Qt, which can be compared to resummation carried 
out directly in Qt space. From the structure of the resummed expressions, we also deter- 
mine the form of nonperturbative corrections to the cross section and implement these into 
our analysis. We obtain a good description of the transverse momentum distribution of Z 
bosons produced at the Tevatron collider. 
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1. Introduction 



The hadronic annihilation cross sections for electroweak boson production (7*, W, Z, H) provide 
important test cases for resummation techniques in QCD. This paper describes the application of 
a simultaneous, joint resummation of threshold and transverse momentum singularities in these 
cross sections. The application to this well-studied case enables us to compare joint resummation 
to the data for Z production at the Tevatron. We regard it as a first step toward extending joint 
resummation to a wider range in kinematics and reactions, and toward a unified description of 
nonperturbative contributions in hadronic reactions. 

At measured transverse momentum Qt *C Q, with Q the pair or boson mass, annihilation 
cross sections have contributions a^\n 2n ~ 1 (Q T /Q 2 )/Q T at each order of perturbation theory. 
Although highly singular at Qt = 0, these terms nevertheless organize themselves into a function 
that describes the Sudakov suppression of the cross section at small Qt [1, 2, 3, 4]. These 
singularities reflect the recoil of the produced electroweak boson against soft gluon radiation. 
This is transverse momentum resummation, or "Qt" resummation, for short. 

Perturbative corrections in factorized inclusive electroweak annihilation cross sections, 

—0^ = J dz J dxi dxj fi/ A (xi, hf) fj/B(xj, Hf) 5 (z - Q 2 /x iXj s) ujI] s (z, Q, //, fi F ) , (1) 

have an analogous behavior at z — 1. In this expression, the /'s are parton distribution functions 
in hadrons A and B, with uj r ^ s the corresponding partonic hard-scattering functions, \x is the 
renormalization scale and /xp the factorization scale. At z — 1, referred to as partonic threshold, 
ujlj S (z) has terms of the form a™ln 2ri_1 (l — z)/(l — z). These singular corrections reflect the 
lack of phase space for soft radiation when the partons have just enough energy to produce the 
observed final state. Like singularities at Qt = 0, they can be resummed to all orders [5, 6]. 

Unlike transverse momentum resummation, the singular functions of threshold resummation 
do not appear as explicit logarithms in the physical cross section, Eq. (1), since they are in con- 
volution with the parton distribution functions. In practice, transverse momentum resummation 
is usually of greater numerical significance at measured Qt than is threshold resummation at 
integrated Qt, simply because the typical value of z in Eq. (1) is generally far from unity. In 
such cases, threshold resummation shows that the effects of distributions singular at z — 1 are 
under control. It also leads to cross sections that are in general less sensitive to the factorization 
scale [7]. For this reason, and because the dynamical origins of the large corrections in both 
threshold and transverse momentum resummations are in soft gluon emission, it is attractive to 
develop a formalism that encompasses both. The necessary analysis for this combination, which 
we refer to as joint resummation, has been carried out in [8] for parton distributions, and in [9] 
for electroweak annihilation as well as QCD cross sections. The effects of resummation are closely 
bound to momentum conservation. The singular corrections associated with soft gluon emission 
exponentiate in the corresponding spaces of moments, impact parameter for transverse momen- 
tum, and Mellin (or Laplace) moments in energy for threshold resummation. The transforms 
relax momentum and energy conservation, while their inverses reimpose it. In joint resummation, 
both transverse momentum and energy conservation are respected. 

In Ref. [10], a preliminary analysis of prompt photon production was carried out from the point 
of view of joint resummation, with an emphasis on the role of recoil in higher-order terms in the 
relevant hard scattering functions. Relatively substantial effects were found in that case. In this 
context, it is important to return to the better-understood example of electroweak annihilation 
to test the joint formalism. 



2 



As mentioned above, there is a close relation between resummation and nonperturbative 
corrections. Taking into account nonleading logarithms and the running of the coupling, resum- 
mation leads in each case to a perturbative expression in which the scale of the coupling reflects 
the value of the transform variable. Because of the singularity of the perturbative effective cou- 
pling at Aq CD , the resulting expressions are, strictly speaking, undefined. A closer look, however, 
shows that singular contributions appear only at nonleading powers of momentum transfer. This 
is an example of how perturbative resummation can suggest the way nonperturbative dynamics 
is expressed in infrared safe hard scattering functions. In effect, perturbation theory is ambigu- 
ous, and the resolution of its ambiguities is, by definition, nonperturbative [11, 12]. As we shall 
review below, these ambiguities manifest themselves as singularities in the integrand of the in- 
verse transforms for both transverse momentum and threshold resummations. Each scheme for 
dealing with these singularities constitutes a specification of perturbation theory, and implies a 
parameterization of nonperturbative effects. We hope that a joint resummation affords a more 
general approach to this problem. 

For each resummed cross section, be it Qt, threshold or joint, one must first specify how 
to define perturbation theory in the transform space, and then, given a well-defined expression, 
how to invert the transform. The most "conservative" approach is to expand the resummed 
expression to some fixed order in perturbation theory [13]. Since the resummation contains 
information on singular terms at all orders, we may in this way get information beyond the order 
at which a complete calculation has been done. To fixed order, the perturbative expression and its 
transform are unambiguous. This approach is applicable to threshold resummation for inclusive 
cross sections, but is not directly useful for phenomenology at measured Qt ~ 0, where the 
cross section is singular. Another approach for Qt resummation is to use the exponentiation of 
leading logarithms in transverse momentum space, and incorporate the finite order improvements, 
including those due to the running coupling and momentum conservation, by a direct expansion 
[14, 15]. A similar application for threshold resummation was developed in [16]. 

Finally, in both Q T and threshold resummation, we may, as indicated above, redefine the 
resummed perturbation theory in transform space, and invert the resulting transform numerically. 
For threshold resummation this has been done by so-called principal- value [11, 17] and minimal 
[18] prescriptions for perturbation theory. Both exploit the analytic structure of the running 
coupling, and redefine transform integrals to avoid the Landau pole. For Qt resummation, 
the commonly-used approach (the ^-prescription) [4, 19, 20, 21, 22] introduces an infrared scale, 
beyond which the running of the coupling is decoupled from the transform variable. More recently, 
Qiu and Zhang [23] have proposed another method, employing a smooth extrapolation of the 
perturbative resummed cross section in impact parameter space to push the Landau singularity 
to infinity. 

In this paper, we will study the jointly resummed cross section derived in [9] , defined by a next- 
to- leading logarithm approximation to the resummed exponent, and an analytic continuation 
of the 6-space contour following [10]. We stress that we could have made other choices, like 
those mentioned above, within the context of joint resummation. We find the approach of [10], 
reviewed below, attractive for its simplicity. It results in a redefined perturbative series with 
no new dimensional scale beyond Aqcd- Phenomeno logical parameters then appear only as 
nonperturbative power corrections. In this "test case" for joint resummation, we shall find a 
consistency with the data from Z production that we believe is comparable to that of other 
approaches. We hope that this will serve as a starting-point for further development. 

We note that different treatments of the perturbative exponent may differ markedly in pure 
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threshold resummation, due to different treatments of nonleading logarithms [17, 18, 24]. This 
important question applies to joint resummation as well, but we shall not explore these differences 
in this paper. This is because the vector boson cross sections discussed here are in the kinematic 
region described above, where z is far from unity in Eq. (1), and where we expect the effect of 
pure threshold resummation to be small in any scheme. 

We begin our technical discussion with a review of the joint resummation formalism of Ref. 
[9], applied to Z production, including a comparison to the exact 0(a s ) cross section. We 
develop the phenomeno logical evaluation of the resulting expressions in Sec. 3. This includes the 
specification of resummed perturbation theory through contour integrals in the two transform 
spaces, the matching to finite order calculations, and the parameterization of nonperturbative 
effects. In Sec. 4, we compare the resulting expressions to CDF and DO data, fitting the necessary 
nonperturbative parameter. Following this, we give a few preliminary conclusions on the status 
and prospects of the joint resummation program. In Appendix A, we present some well-known 
formulas that are relevant for our calculation, and in Appendix B, we sketch the simple numerical 
method that enables us to evaluate resummed moment space integrals for parton distribution 
functions that are specified only in x-space. 



2. The jointly resummed cross section 

Within the formalism of [9], the jointly resummed cross section da rcs / dQ 2 dQ\ for electroweak 
annihilation is obtained as a double inverse Mellin and Fourier (impact parameter) transform: 

d( 7 T B = V (0) f <^i T -N f d% jQ T -b 

V a Jc N 2ni J (2nY ' 



dQ 2 dQ\ V a Jc N 2m J (2tt) 5 

x C a/A (Q,b,N,n,n F ) exp[E^(N,b,Q,fi)] C a/B (Q,b,N,fj,,fi F ) . (2) 

Here, r = Q 2 /S, and is a normalization containing the appropriate electroweak charges 
occurring in the basic underlying process ad — » V . For completeness, and because this is the 
process we will study numerically in this paper, we give for Z production in Appendix A. 

Compared to the expression given in [9], we have brought the cross section into a form 
that is closer to the standard one [4] in Q T resummation. Resummation of soft-gluon effects is 
achieved through the flavor-diagonal Sudakov exponent E^(N, b,Q, //), while the "C-coefficients" 
contain the parton distribution functions and provide resummation of additional logarithms of 
soft-collinear and collinear origin. We will discuss the various terms in turn. 



2.1. The Sudakov exponent 

We now develop an expression for the exponent E^?, valid to next-to-leading logarithms (NLL) in 
N and b, based on the results of Ref. [9] . Compared to that reference, we will absorb some terms 
in the exponent that are associated with parton evolution into the C-coefficients. We begin with 
the exponent derived in [9] for the eikonal approximation to aa annihilation to an electroweak 
boson, E^. It is given by: 



E^(N,b,Q,fi,fi F ) 



rQ- 

7„ 



A a {a s {k T )) 



Ubk T )K { 



(2Nk 7 



+ In 



'Nk 



T 



4 



2 In (TV) /f ^fA a (a s (k T )) . 



Here, Jo and K are the usual Bessel functions, and we define 

N = iVe 7E , 

with the Euler constant. The function A a (a s ) is a series in a s , 



7T V 7T / 



(3) 



(4) 



(5) 



in terms the familiar coefficients [25]: 



= C a (C q = C F , C g = C A ) 



A?) 



a 



! ' 1 18 6 / 9 



(6) 



Dependence on the renormalization scale is implicit in Eq. (3) through the expansion of a s (/cr) 
in powers of a s (fi). 

Following Refs. [9, 4, 26], we approximate the exponent in Eq. (3) by a "minimal" form that 
is accurate to next-to-leading logarithm in both transform variables: 



E^(N,b,Q,fi^ F ) 



Q 2 /x 2 K 



r A a (a s (k T )) In ' 



Q 



-2 In (N) ^fA a (a s (k T )) . 



Here, the function x(iV ', b) organizes the logarithms of N and b in joint resummation, 

1 + r]b/N 

where r] is a constant, and where, by analogy to N of Eq. (4), we have defined 

b = bQe^/2 . 



(7) 



(8) 



(9) 



To anticipate, we will choose rj — 1/4 below. We will discuss the particular expression (8) for x 
in the next subsection. 

We now regroup Eq. (7) as 



Ef*(N,b,Q,n,H F ) 



Q dk 2 T A , f -i w -\ (k T 



2 /x 2 14 



A a (a s (k T )) In 



Q 



- 2 In (N 



Q 2 /x 2 dkl 



n% k\ 



A a (a s {k T )) 



(10) 



Next, we make contact with standard Qt resummation by writing 



H" / 

J U 



Q 2 dk 



T 



Vx 2 K 

Q 2 /x 2 dj4 



A a (a s (k T )) In 



B a (a s (k T )) 



- 2A a (a s (k T )) In (iV) - B a (a s (k T )) 
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Here we have introduced 



7T \ 7T / 



with 



BP = 



Bf = ~{l\C A -AT R N F ) . 



(12) 



(13) 



Eq. (11) follows from Eq. (10) to NLL accuracy in N and b. The coefficients B^f> are also 
known [19, 27], but contribute only beyond NLL. 
The first term in Eq. (11), 



QVx 2 k 2 T 



A a (a s (k T )) \ n [^-)+B a (a s (k T )) 



(14) 



has the classic form of the Sudakov exponent in electroweak annihilation, the only new ingredient 
being the quantity x that depends on N and b and represents the joint resummation. As shown in 
Sec. 3.1 of Ref. [9], the term with B a accounts for the difference between the eikonal approximation 
and the full partonic cross sections in the threshold region. Eq. (14) will be our choice for the 
exponent E^? in Eq. (2). Its expansion in leading and next-to- leading logarithms gives 



hM((3) + 2hW( ( 3,Q,ri, 



where 



hP(J3,Q,n) 



2nb 2 
2nbl 



(15) 



(16) 



2 v fJ 1-2(3 



+ 



2tt6 



ln(l - 2(3) 



+ 



2nb 



2(3 



1-2(3 



+ ln(l - 2(3) 



(17) 



In these equations, 



A 

P 



b Q a s (n) In (iV) , 
b a s (ii) ln(x) , 



and 



bo = 



11C A - AT R N F 
12tt 



bi = 



17C\ - 10C A T R N F - 6C F T R N F 
24tt 2 



(18) 



(19) 



In order to interpret the second term in Eq. (11), we note that the combination in its square 
brackets, 

-2 A a (a s ) In (N) - B a (a s ) , (20) 

corresponds to the leading (logarithmic plus constant) terms at large N in the one-loop diag- 
onal q — > q (or, g — > g) anomalous dimension. To be more precise, the way we have obtained 
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Eq. (11), the second term in Eq. (11) matches the anomalous dimension to NLL accuracy only, 
since the constant (in N) part of the two-loop anomalous dimension is not identical to the cus- 
tomary [19, 20, 21, 22] coefficient (even though it is related to it [27, 28]). As mentioned 
above, contributions related to enter only beyond NLL and are outside the presently devel- 
oped reach of the joint resummation formalism. We are therefore indeed free to associate the 
combination in Eq. (20) with the leading terms in the anomalous dimension. It then follows that 
the second term in Eq. (11) represents the evolution of the parton densities from scale jUj? to 
scale Q/x in the large- N limit, that is, near threshold. Note that indeed all dependence of the 
exponent on /ip is contained in this term. The evolution term we have identified in Eq. (11) 
will become part of the functions C a /n introduced in Eq. (2). We shall therefore pursue it further 
when discussing the C a /H in subsection 2.3. below. 



2.2. The function X (N, b) 



We have defined x{N, b) in Eq. (8). There is an element of choice in the actual form of x(N, b), the 
only requirement being that the leading and next-to-leading logarithms of N and b are correctly 
reproduced in the limits N — > oo or b — > oo, respectively. In ref. [9], the somewhat simpler choice 



x(Ref. [9]) = b + N 



(21) 



was made. While this is a legitimate option, we found it to be less convenient for phenomenolog- 
ical studies. The reason for this is that this form of \ introduces sizable subleading terms into 
perturbative expansions of the resummed exponent, which are not present in full fixed-order per- 
turbative results. For instance, expanding the exponent E^?(N,b,Q, /i) in Eq. (15) to 0(a s ([/)) 
one finds: 

E aa (N,b,Q,H) 



cxp 



TT 



ln (x) - o ln (x) 



(22) 



If we are assuming that b ^> N, relevant at small Qt far away from threshold, then this is 
approximately 



cxp 



Eg(N,b,Q,ii) 



1 _*«m Cf 



TT 



3 N - 

In 2 (b) - -ln(6) +2^1n(6) + 



(23) 



On the other hand, the fixed order 0(a s ) result for the partonic cross section for flavor a is given 
(atQ r ^0)by [2]: 



dQ 2 dQ 2 T 



a 



(0) 



a s (jd) 



C h 



TV 



HQ 2 /Q\ 

Qt 



+ O (HQVQt)) 



(24) 



where we have indicated the functional form of the first Qr-suppressed correction. It is easy 
to show that upon Fourier transformation of the 0(a s ) in Eq. (23) back to Qt space, the first 
two terms ln 2 (6) — §ln(6) reproduce the first two contributions to da/dQ 2 dQ\ in Eq. (24). The 
term oc ln(fo)/fo in Eq. (23), however, yields a subleading contribution to the cross section, which 
is of the form \n(Q^/Q 2 ) /Qt, that is, down with respect to the leading and next-to-leading 
logarithms, but more singular than the first suppressed correction to the fixed order cross section 
in Eq. (24) which is just oc lnlQ^/Q 2 ). In other words, the choice (21) introduces new dependence 
of the resummed cross section on Qt, not present in the cross section calculated at fixed order. 
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Even though this affects only subleading, integrable terms, which are beyond the reach of our 
resummation anyway, this mismatch between Eqs. (22) and (24) at 0(1/ Qt) produces a spurious 
logarithmic singularity in da/dQ 2 dQ?, which is the cross section we will compare to the data. 
This problem is avoided to all orders by choosing \ as hi Eq. (8) with any rj > 0. 

The function x hi Eq. (8) is only a slight modification of Eq. (21), but it has the property that 
at large b corrections to the leading term are suppressed as 1/b 2 . In this way, no dependence of 
the form \n(Q T /Q 2 )/Q T can arise in the cross section in Q T space. Obviously, the general limits 
for b — > oo and N — > oo are the same as in (21). It also turns out that this form of x(N, b) leads 
to analytic properties of the exponent in Eq. (15) that are consistent with the method described 
below for performing the inverse transform in b. 

We will analyze expansions of our final resummed cross section to one loop in more detail in 
subsec. 2.4.. Before doing so, we need to specify the coefficients C a /H(Q,b, N, /i, fip) of Eq. (2). 

2.3. The coefficients C a / H (Q, b, N, /i, /ip) 

The coefficients C a /n(Q, b, N, //, /ip) hi Eq. (2) are chosen to correspond to the jointly resummed 
cross section in [10] for large N and arbitrary b, and to Qt resummation for b — > oo, N fixed: 

C a/H (Q, b, N, /i, fi F ) = J2 C a/J (N, a s (fi)) E jk (N, Q/ x , fi F ) h/ H {N, pt F ) . (25) 

Here the fj/n{N, ftp) are again the parton distribution functions for hadron H at factorization 
scale hf- In principle, by analogy with standard Qt resummation [4], the scale for the strong 
coupling in the C a /j (N, a) would be Q/x] however, at the NLL level we are considering here, it 
is legitimate to choose the "large" renormalization scale // ~ Q and to expand the C a u (N, a) to 
a finite order in a s . To first order, matching Eq. (25) to the large-6 behavior of the Qt- resummed 
cross section [2, 4, 19, 28], one has 

C q/q (N, a.) = l + gcv(-8 + 7r 2 + N J +1) ^j = C m {N,a a ) , (26) 
a 1 

c ^ N - a ' ] = t WTTWTJ) = c * (iv ' a ' ) ' (27) 

Note that the coefficient C q / g is off-diagonal in flavor. It is indeed a well-known feature in Q T 
resummation [4, 19, 2, 28] that such non-diagonal terms also contribute to singular behavior at 
Qt = 0. On the other hand, they do not incorporate singularities at threshold, which is visible 
from the fact that C q / g is suppressed at large N . 

The evolution matrix £ (N, Q/x^f) in Eq. (25) results from the second term in Eq. (11) that 
we chose to absorb into the C coefficients. Compared to the large- N limit used in that equation 
and relevant near threshold, we can make an improvement here and replace the leading- iV part 
of the diagonal anomalous dimension in Eqs. (20) and (11) by the full anomalous dimension 
relevant for the scale evolution of parton densities, containing also all terms subleading in iV: 

-2A a (a s ) In (n) - B a (a s ) — > lN {a s ) . (28) 

At NLL level, where contribute, we need the first two terms in the perturbative 

expansion of the anomalous dimension [29], 7^ = ^7yv' ) + (^f) 1$ ■ Since the 7^ are (in 
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general) matrices, this procedure introduces terms that are parton non-diagonal, and thus leads 
to the matrix structure of £ in terms of an ordered exponential. As mentioned above, the 
interpretation of £ is simply the evolution of the parton densities from scale \ip to scale Q/x, 
within the NLL approximation. In this way, it leads to a resummation of collinear logarithms, 
some of which are associated with partonic threshold, that is, are also proportional to IniV, 
while others are suppressed by 0(1/ N) or more, due to partonic mixing. Such a procedure 
is also familiar from standard Qt resummation [4, 28], and the substitution in Eq. (28) thus 
provides a natural extension of our formalism for joint resummation away from threshold. We 
may also interpret C a /H{Q,b, N, /i, fip) as the Fourier-Mellin transform of a generalized parton 
distribution, at measured transverse momentum and energy fraction [9]. Eq. (25) then has the 
interpretation of a refactorization of C a /n into a coefficient function C a /j and light cone parton 
distribution functions, fj/H, at the scale Q/x- 

Explicit expressions for the solution of the standard evolution equations for parton densities 
between scales [ip and Q/x can be found in [30, 31] and can be used to determine the elements 
of the matrix £. To achieve the exponentiation of the evolution terms - despite the fact that the 
matrices 7^ and 7^ do not commute - the iterative procedure derived in Ref. [31] is particularly 
useful. The parameter that governs the evolution between scales hf and Q/x is 

where the right-hand-side is the expansion to NLL accuracy, consistent with our approximation. 
It should be emphasized that in the above expression the scale hf appears through a single explicit 
logarithm that will serve to approximately cancel the [ip dependence of the parton distributions 
in Eq. (25), resulting in a decrease in /z^-dependence for the final resummed cross section in 
Eq. (2), as compared to a fixed-order calculation. 

In all of these procedures, working in Mellin-iV moment space is a great convenience, be- 
cause it enables us to explicitly express the evolution between the scales fiF and Q/x in terms 
of the parameter \n(a s (iiF) / ®s(Q / x)) in Eq. (29). In this way, we avoid the problem normally 
faced in Q T resummation that one needs to call the parton densities at scales far below their 
range of validity, so that some sort of "freezing" (or related prescription) for handling the par- 
ton distributions is required. As is evident from Eq. (25), we only need the parton distribu- 
tion functions at the "large" scale hf ~ Q, whereas normally in Qt resummation the product 
J2k£jk ^N,Q/b, jj,F^j fk/n(N, hf) is identified with fkj}j(N,Q/b). An organization of the C co- 
efficients in a form similar to ours was first proposed in [28]. We finally note that the moment 
variable N and, as will be discussed below, also the impact parameter b are in general complex- 
valued in our approach, so that it is even more desirable to separate the complex scale Q/x 
from that in the parton densities. In this way, it is not even necessary (albeit convenient) to 
have the parton densities in Mellin-iV moment space, as provided in the code of [32]. In fact, we 
can generalize our analysis to any set of distribution functions, even if specified only in x space. 
Details are discussed in Appendix B. 

2.4. Finite-order 0(a s ) expansions of the resummed cross section 

In this section, we compare expansions of our jointly resummed cross section to "exact" ex- 
pressions for the electroweak annihilation cross section at 0(a s ). We discuss first the limits 
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N — > oo, 6 = and 6 — > oo, A" fixed separately. In these cases, all results can be given 
analytically. 

The limit iV — > oo, 6 = corresponds to pure threshold resummation for the total cross 
section. It is realized by integrating over Qt in Eq. (2), which sets 6 = there. Expansion 
of Eq. (2) to 0(a s ), using Eqs. (15), (26), (27) and the parton distribution functions evolved 
according to Eq. (29), then gives for the partonic cross sections in the qq and qg scattering 
channels: 



r (0) 



a a' 77" °f { -4 ln z N + 6 In JV - 8 + n z + 
2n I 



+ 



N(N + 1) 



AT(AT + l) 
+ 3-4S' 1 (JV) 



r (0) «» T 



N 2 + N + 2 
q ' 2n ±R \ N(N + 1)(N + 2) 



-2 IniV + In 



A, 



2 In N + In ^ 

9 



(iV + l)(iV + 2) 



(30) 



where Si(N) = J2f=i j 1 = 4>(N + 1) +7#, with ■?/> the digamma function. Since we are interested 
in the near-threshold region, we can expand this further to the large- iV limit. Using ip(N + l) = 
In N + 1/(2N) + 0(l/iV 2 ), we find: 



r 19 



2tt 



C F <U In 2 + 4 



IniV 

~aT 



- 8 + 7T 2 + 



^ 2^iV 



-2 In N + In [ %- 



+ 



3 41niV 

N 

(\nN\ 



In 



A, 



+ o 



In AT 

Iv 2 " 



A" 2 



(31) 



The Mellin moments of the "exact" 0(a s ) partonic cross sections can be found in [33] and read: 



w exact 



2tt 



N(N + 1 



N 2 (N+l) 2 



-8 + 



4tt 2 



+ 



N(N + 1) 



+ 3-4Si(AT) 



In 



^ exact 



2 A^ 2 + A^ + 2 s ^ + + 11N 3 + 22N 2 + UN + 4 



N(N + 1)(N + 2) 



N 2 (N + 1) 2 (N + 2) 2 



+ 



N 2 + N + 2 



N(N + 1)(N + 2) 
At large N, this gives 



In 



u exact 



cr(°)^C F <|41n 2 iV + 4- 



qg 

exact 



2tt 

<x (0) — Tr — 
9 2tt TV 



A" 



-21niV + In 



3 4 In A" 

N 

'lnA>\ 



{i*F, 



(32) 



+ o 



InN 



N 2 



N 2 



(33) 



Comparing Eqs. (31) and (33), we see that at large A" the expansion of the resummed cross 
section correctly reproduces the 0(a s ) result, including even all terms that are down by 1/N. 
The only difference between (31) and (33) is in the term oc 7r 2 . This difference is due to our 
choice of the coefficient C q / q in Eq. (26), for which we employed a form that is more standard 
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in Qt, rather than in threshold, resummation. A closer inspection of our resummed eikonal 
exponent, Eq. (3), reveals that its Bessel functions result in different contributions oc n 2 in the 
two limits N — > oo and b — > oo, just as needed to explain the deficiency between Eqs. (31) and 
(33). In other words, we could modify our expansion of the exponent into logarithms somewhat 
(by suitably redefining x), so that Eq. (31) would automatically have the correct coefficient of 
7T 2 . On the other hand, the terms associated with ir 2 are beyond NLL, which is the scope of the 
present anlysis, and we therefore do not implement this change here. 

It is worth pointing out that the reason why we correctly reproduce all terms suppressed as 
1/N in Eqs. (31) and (33) is our treatment of evolution in Eq. (28). As was discussed in [34], 
the leading \n(N)/N terms are associated with collinear non-soft emission; it is therefore natural 
that they can be generated from evolution of the parton densities between the scales fip and 
Q/N, as embodied in Eq. (29) at b = 0. In this way, our joint resummation corectly includes 
the leading a k s ln 2fc_1 N/N terms to all orders. Because threshold resummation has a relatively 
modest effect for vector boson production in the kinematic region explored at the Tevatron, we 
leave for future work a more complete comparison of our resummed expression to fixed order, 
beyond these 0(a s ) considerations. 

In the limit b — > oo, N fixed, our formulas smoothly turn into those for standard Qt resum- 
mation. For the one- loop expansion of the jointly resummed cross section we find 



• r qq 



,-(o) 



a s 
27 



C F <^-41n 2 6 + 61n6-8 + 7T 2 + 



N(N + 1) 



(34) 



+ 



+ 3 



a qg 



_N(N + 1) 

r n 2 + n + 2 

2^ lR \N(N + l)(N + 2) 



r (o) ^ T 




21n6 + ln 



+ 



+ 0{\nb/b 2 ) 



(N + l)(N + 2) 



+ 0(}nb/b 2 ) , 



in full agreement with the expressions for the large-6 limit of the "exact" 0(a s ) result derived 
in [2]. 

We do not present closed expressions for arbitrary large N and b, but we can easily compare 
numerically the exact 0(a s ) result with the expansion of Eq. (2) to 0(a s ). Fig. 1(a) shows the 
fractional deviation 



A = 



rf(7 fixcd(l) d a exp(l) 



/ 



rf(T fixed(l) 

dQ T 



(35) 



dQ T dQ T 

where da &xcd ^ /dQ T is the "exact" 0(a s ) cross section and da cxp ^ / dQr denotes the one-loop 
expansion of the resummed expression. Fig. 1(b) compares da &xcd ^/dQ T and da cxp ^/dQ T 
individually. Note the excellent agreement in the region Qt < 10 GeV where resummation is 
necessary. Beyond 10 GeV, the agreement is naturally less exact but still good. This is the region 
where matching to finite order is appropriate, to which we will turn now. 



3. Inverse transforms and matching 

In this section we take the remaining steps necessary to apply the joint resummation formalism 
phenomenologically. This includes specifying a prescription for performing the inverse integrals 
in (2), as well as implementing a procedure for matching resummed and finite-order results. In 
addition, we will consider nonperturbative effects resulting from the strong coupling at small 
momentum scales. 
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Figure 1: (a) Fractional deviation A (as defined in Eq. (35)) between the "exact" 0(a s ) result and 
the 0(a s ) expansion of the jointly resummed cross section. We consider here Z boson production 
at the Tevatron; the cross section has been integrated over 66 < Q < 1 16 GeV. (b) Comparison 
of da^ ed ^/dQ T (solid) and da cx ^/dQ T (dashed) on absolute scale. 

3.1. Inverse transforms 

When performing the inverse N and b transforms, special attention has to be paid to the singu- 
larity in the resummed exponent, Eqs. (15)— (17), at (5 — 1/2. For rj = 1/4, the singularity occurs 
for 

- - - N 

X (N, b) = b+ - x = exp [l/(26 a 8 ( A i))] = p L , (36) 

AN 

and is a manifestation of the Landau pole in the strong coupling. The exponent is also ill- 
defined when x = and infinity, i.e., at b — —2N and b = —AN, respectively. We note that 
the choice 77 = 1/4 is made simply to reduce x = to a linear relation between b and N. This 
is not an essential simplification, but it makes the following analysis slightly more convenient. 
In the following we sketch the application of the method of [10] to this jointly resummed cross 
section. As far as the Landau pole is concerned, it is simultaneously in the spirit of both the 
"minimal" prescription proposed for pure threshold resummation in [18], and the "principal value 
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resummation" , described in [11]. 

The contour for the inverse Mellin transform is chosen to be bent at an angle with respect 
to the real axis and is parameterized as follows (see Fig. 2): 

N = C + ze^ , (37) 

where the upper (lower) sign applies to the upper (lower) branch of the contour, with < z < oo 
(oo > z > 0). For (j) > 7r/2, this results in an exponentially convergent integral over N in the 
inverse transform, Eq. (2) for all r < 1 [16, 18]. For the moment, we do not specify the constant 
C, except that it has to lie to the right of the rightmost singularity of the parton distribution 
functions. At any finite order in perturbation theory, all values of C > 0, and all it > > 7r/2 
are equivalent. In the resummed cross section, however, the singularity at x — Pl introduces a 
power-suppressed ambiguity in the transform, which we resolve by choosing C < pl [18] . 




As mentioned above, the position of the Landau pole depends on both and b. The d 2 b 
integral in Eq. (2) can be written as 

d 2 be iQT - b f(b) = 2tt ^ dbb J (bQ T ) f(b) , (38) 

where Jo is the Bessel function. Hence, as it stands, the integration over b lies on the positive 
real axis. This would imply that one would never be able to avoid hitting the Landau pole, 
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since for any choice of the parameter C in Eq. (37) there is a b > for which x(N,b) = Pl- 
This problem is of course well-known from standard Qt resummation. The procedure usually 
adopted is to prevent b from becoming too large by evaluating the resummed cross section at 



fe* = b/yj 1 + fr 2 /& max [4, 3, 19, 20, 21, 22], at the expense of introducing a new parameter 6 max . 
To avoid introducing a new parameter, we treat the b integral in a manner analogous to the 
N integral above, avoiding the Landau singularity on a contour that produces an exponentially 
convergent integral for all Qt > 0. 

Were the Landau pole not present we could, instead of performing the b integral along the 
real axis, use Cauchy's theorem and divert it into complex b space along either of the two solid 
lines in Fig. 3, under the condition that the integrand falls off sufficiently fast at large |6| and 
that there be no contribution to the integral at infinitely large real part. To achieve this, we have 
to split up Eq. (38) as [10] 

rOO POO 

2tt dbb J (bQ T )f(b)=n dbb[h 1 (bQ T ,v) + h 2 (bQ T ,v)} f(b), (39) 

JO JO 

where we introduce two auxiliary functions hi i2 (z,v), related to Hankel functions and defined in 
terms of an arbitrary real, positive parameter v by integrals in the complex #-plane [35]: 

1 [■—■K+iVTT 

hi(z,v) = / d9e~ izshie , 



TT 



IVTV 



1 r—ivTT 

h 2 (z,v) ee — / d9e~ izshld . (40) 

TT J IT+iviT 

For hi, we parameterize 9 = —ivTT + xoTT(—l + 2iv) (0 < xq < 1), while for h 2 , 6 = — iv tt + X0Tt(1 + 
2iv) (1 > xe > 0). The hi i2 become the usual Hankel functions Hi j2 (z) in the limit v — > oo. 
We note that this convergence to the Hankel functions is extremely rapid, since the dependence 
on the variable v is suppressed by the exponential of an exponential for all finite z. The 
are finite for any finite values of z and v. Their sum is always: hi(z,v) + h 2 (z,v) = 2J (z), 
independent of v. The utility of the h- functions is that they distinguish positive and negative 
phases in Eq. (39), making it possible to treat the b integral as the sum of the two contours in 
Fig. 3, the one associated with hi (h 2 ) corresponding to closing the contour in the upper (lower) 
half plane. 

The virtue of this technique for the b integration is that we can choose the contours to avoid 
the Landau pole. We simply need to make sure that the b contour never intersects the trajectories 
defined by x(iV, b) = p L , shown by the two light solid curves in Fig. 3. As mentioned earlier, 
singularities also arise for b = —2N and b = —4N. These contours are shown in Fig. 3 by the 
dotted line and the dash-dotted line, respectively. Parameterizing the upper b contour as 

Cl - b ~\ bc-te-** (0<t<oo) (41) 

and the lower one as 

r ■ b - { 1 (o<t<b c ) 

° 2 - b -\ b e -U?*> (0<t<oo) ' (42) 

we choose the parameters b c and <ftb such that none of the branches intersects any of the "forbid- 
den" lines in Fig. 3. A typical choice is also shown in Fig. 3 by the thick solid lines. The parame- 
ters C in Eq. (37) and b c in Eqs. (41), (42) are arbitrary as long as < ( 1+b ^ Q / 8 c + ^§~) e7E < Pl- 
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Figure 3: Choice of contour for b integration (thick solid lines). The straight sections of the 
contour from to b c are to be interpreted as on the positive real axis. The remaining curves 
represent lines of singularity discussed in the text. 



In this way, our full expression for the cross section in terms of inverse transforms, Eq. (2), be- 
comes 



dQ 2 dQr „ ° JCn \ J a 47r 

dbb 



+ j —h 2 (bQ T ,v)wii(Q,b,N,n,LiF) , (43) 



with 



wj&(Q,b,N,»,» F ) = exp[E™(N,b,Q,»)' 

x S C a/j (N, a 3 (n)) £ jk (N, Q/ X , Hf) fk/A(N, (i F ) 

x E Qtf W «<M) % (N, Q/x, M h /B (N, im f ) , (44) 

jk 

for initial hadrons A and -B. 

This choice of contours in complex transform space is completely equivalent to the original 
form, Eq. (38) when the exponent is evaluated to finite order in perturbation theory. It is a natural 
extension of the iV-space contour redefinition above [16, 18], using a generalized "minimal" [18] 
exponent, Eq. (15). As we stressed earlier, joint resummation with its contour integration method 
provides an alternative to the standard b space resummation. Joint resummation has built-in 
perturbative treatment of large b values, eliminating the need for a 6* or other prescription for 
the exponent, or for a freezing of the scale of parton distributions at large b or low Q T . 

To examine the relevance of the large-6 contributions, we now compare the jointly resummed 
Qt distribution matched to the 0(a s ) perturbative result [2] and the Qt space resummed [15] 
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Figure 4: Q T distribution for Z production at y/s =1.8 TeV calculated using the narrow width 
approximation (Q = Mz). The Qt space method result and the joint resummation method 
result are matched to the cross section at 0(a s ). The CTEQ5M [36] parton distributions have 
been used. 



distribution also matched to the 0(a s ) distribution (see the next section for details of the match- 
ing procedure). The Qt space resummation formalism originates from b space resummation and 
can be viewed as a very good approximation of the latter. Unlike the standard b space technique, 
direct Qt resummation, like the contour method just described, yields a result even without any 
nonperturbative input for nonzero values of Qt- The two approaches are compared in Fig. 4 
for the case of Z boson production at Tevatron energy. As expected, the two distributions differ 
mostly at the very low-Qr end of the spectrum. 



3.2. Matching to finite order 

A resummed Qt distribution provides theoretical predictions for the small Qt region; in the 
large Qt regime one mainly relies on fixed-order perturbation theory. Thus a description of the 
intermediate Qt region requires a consistent matching between the two results that avoids double 
counting. For the joint resummation we adopt a matching prescription first proposed in standard 
b space resummation [20]: 

da da res da cxp M da &xcd ^ 

dQ 2 dQ T ~ dQ 2 dQ T ~ dQ 2 dQ T + dQ 2 dQ T ' ^ ^ 
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where da rcs / 'dQ 2 dQ\ is given in Eq. (2) and, as before, da cxp ^ /dQ 2 dQ^ denotes the terms 
resulting from the expansion of the resummed expression in powers of a s (/x) up to the order k 
at which the fixed-order cross section da 5xcd ^ /dQ 2 dQ\ is taken (in practice, k — 1 (see [2]) or 
k = 2 (see [37])). 

Alternatively, one can formulate the matching procedure in the following way [4, 2]: 

da da ies 

— W = 7-02,7/02 + ^finite , (46) 



dQ 2 dQ 2 T dQ 2 dQ\ 

^finite standing for the finite part of the fixed-order distribution, i.e. the part remaining after 
the singular behavior oc a s \n m (Q 2 /Ql)/Q^, (2k - 1 > m > 0) at fcth order of perturbation 
theory (in our case, k — 1,2) has been taken out of the fixed-order cross section. For pure Qt 
resummation, there is no difference between the two ways of matching in Eqs. (45), (46), as long 
as G?<7 cxp( W / dQ 2 dQ\ coincides identically with the singular part of the fixed-order cross section. 
For our joint resummation, however, rfa cxp ^ / dQ 2 dQ 2 T also contains terms that are non-singular 
in Qt (as in Eq. (23)), which implies that the matching is most naturally performed according 
to Eq. (45) in conjugate (6, N) space, in order to avoid double counting. 

Notably, joint resummation with the choice of x{N, b) in Eq. (8) and the matching prescrip- 
tion (45) returns a positive cross section even in the high Q T regime. This is in contrast to pure 
Qt resummation, which notoriously yields a negative answer at large Qt- To avoid the latter 
feature, usually an additional switch between the matched cross section and the fixed-order result 
is implemented in pure Qt resummation. A nice feature of joint resummation is that this is not 
necessary here; Eq. (45) is all we need all the way to Qt of order Q. (Of course, for Qt Q 
a further resummation is necessary [38].) The improved large-Qr behavior of the matched cross 
section can be traced back to the behavior of the jointly resummed part da res / 'dQf ' J dQ T itself, 
which also remains positive for all values of Qt- This in turn results from the behavior of the 
function x{N,b) at the small values of b relevant to high Qt- 

A simple argument shows how the small-6 behavior of x ensures that the cross section remains 
positive at large Qt for the jointly resummed cross section, and also why it goes negative in pure 
Qt resummation in impact parameter space. Consider first pure Qt resummation, which is 
equivalent to the choice x — b- For any the one-loop expansion of the resummed exponent 
in impact parameter space is given by Eq. (22). When x — b simply (N = in Eq. (23)), this 
expression is just as singular at b = as for b — > oo [39]. This produces a spurious suppression for 
b small compared to 1/Q, which is present as well in the full exponent. Because the b integral is 
dominated by b ~ 1/Qt, this effect is unimportant when Q T is small. On the other hand, when 
Qt is large, of the same order as Q, the suppression of small b removes an important contribution 
from the b integral. Along the deformed contours of Fig. 3, the situation is particularly simple. 
Referring to Eq. (43), the sum of the functions hi(bQT, v) and hi^bQT, v) starts out at 2 for b = 0, 
but both begin to oscillate and to decrease exponentially, when b ~ 1/Qt- A spurious suppression 
times W^(Q,b,N, H,Hf) for b < 1/Q ~ 1/Qt will thus eliminate positive contributions to the 
transform, leaving over negative contributions, in the limited range of b in which hi and h 2 are 
themselves of order unity. If, on the other hand, we pick x as in Eq. (8), then we ensure that the 
spurious suppression at small b characteristic of pure 6-space resummation is absent. The result 
is a positive cross section for large Qt- 
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3.3. Nonperturbative input 



In most applications of b space resummation, the perturbative component is supplemented by a Q- 
dependent Gaussian in impact parameter space, to tune the overall influence of nonperturbative 
dynamics [4, 19, 20, 21, 22]. Joint resummation implies similar effects in the Q T distribution 
at small Qt- The starting point is the full NLL exponent [9] given in Eq. (3). We expand the 
Bessel functions in Eq. (3) to derive an explicit form for the leading b- and iV-dependence of 
the resummed exponent, times an integral of the anomalous dimension A a over soft transverse 
momenta, defined by a cutoff A > Aq CD , 

/ fr2 27V 2 \ r x2 / O \ 

E£(N, b, Q, /x, » F , A) ~ (-- + —j jf dk 2 T A a (a s (k T )) In ( JjL-J . (47) 

In this way we derive a standard Gaussian form in b for a multiplicative, nonperturbative smearing 
function in b space, to account for effects from small kr- The new feature here is that the 
coefficient of b 2 is essentially identical to that of the threshold-related power correction (N/Q) 2 . 
In our numerical applications below, we will be rather far away from threshold, and as a result 
we will be mainly sensitive to moderately small N. Hence we will retain only the term ~ b 2 
in Eq. (47). Thus our nonperturbative smearing function reverts to the purely Gaussian one 
usually used in Q T resummation, and we will make the following replacement in the resummed 
cross section: 

E P J(N, 6, Q, pi) — E P J(N, 6, Q, pi) - gb 2 , (48) 

where g is a parameter to be determined by comparison to data. Note that Eq. (47) implies that, 
as usual, g has a component that depends logarithmically on Q [4, 19, 20, 21, 22, 40]. 

4. Vector boson production in the framework 
of joint resummation 

With the developments discussed above, the joint resummation formalism becomes a practical 
tool for the description of electroweak annihilation. Here we consider Z boson production at the 
Tevatron collider. Recent data on the Qt distribution of the produced Z bosons are available 
from both CDF [41] and DO [42] experiments and have reached a good level of precision. The 
overall normalization of the data will be treated as a free parameter in our analysis, and will be 
varied within the quoted experimental errors. 

All results are obtained using CTEQ5M parton distribution functions [36], in the manner 
described in Appendix B. We choose the factorization and renormalization scales pi = pip = Q. 
The numerical values of the electroweak parameters we use are as follows: Mz = 91.187 GeV, 
sin 2 9\y = 0.224, V z = 2.49 GeV. Note that the experimental data sets have been integrated over 
finite regions in Q, 66 < Q < 116 GeV for CDF and 75 < Q < 105 GeV for DO. The inverse 
transforms are performed as described in Section 3.1., with the following choice of the contour 
parameters in Figs. 1 and 2: <p = <p b = 25/32 tt, C = 1.3, b c = 0.2/Q. Of course we are free to 
choose these parameters differently, as long as they are such that the structure of the contours 
as depicted in Fig. 3 is maintained. Finally, as mentioned above, we choose rj = 1/4 in Eq. (8). 
We have checked that the result is quite insensitive to this choice, at the order of a percent when 
r] is changed by a factor of 2 at Q T = 4 GeV. 
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As we pointed out before, an attractive feature of the joint resummation with transforms 
defined as above is that we can obtain predictions that have no dependence on any additional 
nonperturbative parameter. Results of this form are shown by the dashed lines in Figs. 5 and 6. 
In both cases, we have adjusted the normalization of the theory curve so that the x 2 °f the 
comparison between data and theory becomes minimal. For the CDF data, this normalization 
factor is 1.035, for DO it is 0.96. It is evident from Figs. 5 and 6 that our "purely perturbative" 
predictions correctly reproduce the trend of the data over most of the measured region in Qt, but 
peak at too small Qt- This is no surprise of course, since we expect nonperturbative effects to 
play a non-negligible role at low Qt- For comparison, Fig. 5 also displays the fixed-order (0(a s ), 
dotted lines, and 0(a 2 ), dash-dotted lines x ) results for the cross section, with their well-known 
divergent behavior at small Qt- 

Interestingly, the fixed-order 0(a s ) result misses the data also at large Qt, where it remains 
too low even if the normalization is adjusted within the errors quoted in experiment. Joint re- 
summation, with the choice Eq. (8) for x(-^> b) an d the matching procedure described in Eq. (45), 
adds an important contribution to the cross section also here: the difference da res /dQ 2 dQ T — 
da cxp ^ /dQ 2 dQ T in Eq. (45) remains numerically significant also at large Q T and appears to be 
crucial for bringing the theoretical calculation to the data. As can be seen from Fig. 5, the cross 
section at 0(a 2 s ), dcr fixcd(2) jdQ 2 dQ 2 T , is larger than cfo- fixcd(1) jdQ 2 dQ 2 T . As expected, the difference 
da Tes /dQ 2 dQ 2 T - <i<7 cxp(2) jdQ 2 dQ 2 T is smaller than da XQS j dQ 2 dQ 2 T - da cxp{1) jdQ 2 dQ 2 T , so that the 
full cross section in Eq. (45) depends only little on whether matching is performed at 0(a s ) or 
at 0(a 2 ). We take this feature as an indication that our approach for matching, Eq. (45), is 
justified and reasonable even at large Q T ~ Q- 

To improve the Iow-Qt behavior further, we introduce a nonperturbative function as described 
in Eq. (48). Since most of the cross section comes from the region Q ~ Mz, we neglect the mild 
(logarithmic) dependence of the nonperturbative parameter g on Q. We then fit g to the CDF 
and DO data simultaneously, allowing again the normalizations to vary for the two data sets. 
Since g should be determined from the \ow-Qt region, we include only the data points with 
Qt < 50 GeV in this fit. The results of the fit do, however, not depend much on this choice. 
The optimal result is obtained for g = 0.8 GeV 2 and the normalization factors 1.069 (0.975) 
for CDF (DO). The x 2 for the 42 (20) data points from CDF (DO) included in the fit is 31.3 
(23.4). Even better fits would be possible by using a more refined nonperturbative function, as 
suggested by Eq. (47), resulting in extra parameters. As can be seen from the solid lines in Figs. 5 
and 6, with these parameters a very good agreement between the jointly resummed cross section, 
matched with the fixed-order cross section at 0(a s ), and the data is achieved. We note that our 
nonperturbative parameter g = 0.8 GeV 2 is very similar to that determined in Ref. [23], where 
an extrapolation of the perturbation theory result to large b was made. Ref. [23] argued that the 
b integral in pure-Qr resummation is dominated by the saddle point. Our method and theirs 
lead to somewhat smaller values of the nonperturbative parameter than when a 6* prescription 
is used [21]. 

5. Conclusions 

The transverse momentum distribution of the Z is by now a well-studied problem, to which a 
number of successful analyses have been applied [14, 15, 21, 22, 23]. We have come back to 

1 We have used subroutines of the RESBOS package of Ref. [22] in order to calculate the cross section at 0(a 2 s ). 



19 



this topic because we believe that the method of joint resummation offers additional insight 
on the interplay of perturbative and nonperturbative corrections in this, and other hadronic 
reactions. Joint resummation, implemented as above, provides a convenient definition of the 
perturbative cross section at any nonzero Qt, without the introduction of additional dimensional 
scales (beyond Aq CD ) to define either the perturbative resummation or the parton distributions 
at low scales. 

Treated this way, the jointly resummed cross section retains its original perturbative asymp- 
totic expansion order by order. It also suggests the functional form of nonperturbative corrections. 
Because perturbative and nonperturbative components of a QCD cross section are linked at the 
level of power corrections [11, 12], it will be necessary, and useful, to reanalyze the functional 
and phenomenological aspects of nonperturbative corrections in this formalism, relying on the 
available range of data for the Drell-Yan mechanism. Toward the high-energy side, a further 
application to Higgs production [28, 43] at the LHC will also be of interest. In the same spirit, 
we intend as well to return to the application of joint resummation to semi-inclusive processes 
such direct photon production [9, 10]. 



Acknowledgments 

We are grateful to A. Vogt for valuable discussions on Ref. [31], and for providing a version of his 
code for the NLO evolution of parton densities in which the iterative procedures for approximating 
the ordered exponentials are implemented. We are also thankful to C. Balazs, P. Nadolsky and 
CP. Yuan for helpful communications on the use of the RESBOS code [22], and we thank D. de 
Florian, E. Laenen and J.-W. Qiu for useful discussions. The work of G.S. was supported in part 
by the National Science Foundation, grants PHY9722101 and PHY0098527. W.V. is grateful to 
RIKEN, Brookhaven National Laboratory and the U.S. Department of Energy (contract number 
DE-AC02-98CH10886) for providing the facilities essential for the completion of this work. A.K. 
thanks the U.S. Department of Energy (contract number DE-AC02-98CH10886) for support. 



A Some useful formulas 

For Z boson production one has in Eq. (2) the standard tree- level cross sections (see [44]): 



r (0) = l^g2 

9rS 2 a 

2 Q\Q 2 -Ml) 



K = 



+(af + vf)(al + viy Q> 



(Q2 _ M 2 )2 + M | r | 

V2G F Ml 



Ana 

1 1 
aj = -- , vi = -- + 2sin 2 9 w 



a a = T A a , v a = T*-2e a sm 2 9 w . (49) 
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B Using x-space parton distributions 



The above formulas are directly applicable if the parton densities, including their evolution, are 
treated in Mellin moment space. In this context, it is convenient to use the evolution code of [32] 
which is set up in moment space. In practical applications, however, one may prefer to be more 
flexible concerning the choice of parton densities and be able to make direct use of any (x-space) 
parameterization on the market [36, 45]. One way of achieving this was presented in Ref. [18]. 
Here we propose a new simple method. Let us first rewrite Eqs. (2), (25) as 

d&dQ*. = 2^ ic N dNT ~ N ? MN^ F ) f j/H (N,»F) &*{N,Q,Q t ,m f ) , (50) 

which is obtained after performing the d 2 b integration in Eq. (2). The <r^ es are then resummed 
partonic cross sections, differential in Qt- 

The inverse-moment expression in Eq. (50) is of course identical to an x-space convolution of 
the parton densities with the resummed cross section : 

da res „ r°° dz f 1 dy 2 7 ( r 



dQ> dQi - § r i L f ky ' M h [b ^ ) ^ Q ' Qti M ' (51) 

where z = Q 2 /s, the fi(x,fi F ) are the x-space parton densities, and a™ s (z,Q,QT, (i, fi F ) is given 
by the inverse Mellin transform of the moments &lj S (N, Q, Q T , fx, li f ), 

al] s (z,Q,Q T ,^fi F ) = ^~ J dN z~ N &^(N,Q,Q T , li, li f ) . (52) 

As was pointed out in [18], and as is indicated by the upper limit oo in Eq. (51), the function 
<3\f(z, Q, Q T , fi, Hf) defined in the "minimal" prescription is non-vanishing also at z > 1 due to 
the presence of the Landau pole to the right of the N space contour (see Fig. 2), even though 
it rapidly decreases with increasing z. At z > 1, the angle <p of the A^-space contour has to be 
decreased to below it/ 2 to obtain a convergent result. 

The right-hand side of Eq. (51) in principle allows for using a;-space parton distributions. 
However, a problem arises from the fact that the resummed cross section is highly singular [16] 
at z — > 1 (even though regularized in terms of plus-distributions), which makes the convolution 
with the parton densities numerically very tedious [18]. A convenient way of eliminating this 
problem is to trivially rewrite Eq. (50) as 



2~ri Jc N ^ [(N-l) f i/H (N,Li F )\ [(N-l) f j/H (N,Li F )\ {N _ l)2 

(53) 



dQ 2 dQ\ 2mJcs tT L *^v«v",^yj L v- ^^ v ..,^yj {N _ 1)2 

The Mellin-inverse of aj] s (N, Q, Q T , li, fx F )/(N - l) 2 , 

S l3 (z, Q, Q T , » F ) = — dN z ^-^ , (54) 

is now sufficiently well-behaved at large z thanks to the extra suppression by l/(N — l) 2 . For the 
inverse of (N — 1) fi/ii(N, fip) one finds, making use of the fact that the x-space parton densities 
vanish at x — 1: 

^- f dN x~ n (N — 1) fi/ H (N, /ip) = -^-[xf t (x,fu F )} = F(x,hf) ■ (55) 
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We thus arrive at 



which has good numerical behavior. The standard sets [32, 36, 45] of the parton distributions 
allow taking the first derivative numerically. Depending on the large- N behavior of the resummed 
cross section in moment space, alj s (N,Q,QT, fi, (If), it may be necessary to divide by a higher 
power of iV — 1 in Eq. (54), resulting in higher derivatives of the parton distributions in Eq. (55). 
This turns out to be the case for the gluon-gluon initial state in inclusive Higgs production via 
<7<7-HX [46]. 
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Figure 5: CDF data [41] on Z production compared to joint resummation predictions (matched 
to the 0(a s ) result according to Eq. (45)) without nonperturbative smearing (dashed) and with 
Gaussian smearing using the nonperturbative parameter g = 0.8 GeV 2 (solid). The normaliza- 
tions of the curves have been adjusted in order to give an optimal description; see text. The 
dotted and dash-dotted lines show the fixed-order results at 0(a s ) and O(o%), respectively. The 
lower plot makes the large Qt region more visible. 
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Figure 6: Same as Fig. 5, but compared to the DO data [42]. 
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